Theory of charge sensing in quantum-dot structures 
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Charge sensing in quantum-dot structures is studied by an exactly solvable reduced model and 
numerical density-matrix renormalization group methods. Charge sensing is characterized by the 
repeated cycling of the occupation of current-carrying states due to the capacitive coupling to trap 
states which are weakly coupled to the leads. In agreement with recent experiments, it results in a 
variety of characteristic behaviors ranging from asymmetric Coulomb-blockade peaks to sawtooth- 
and dome-like structures. Temperature introduces distinct asymmetric smearing of these features 
and correlations in the conductance provide a fingerprint of charge-sensing behavior. 
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PACS numbers: 73.23.Hk,71.15.Dx,73.23.-b 

Introduction. — Within the orthodox picture of the 
Coulomb blockade regime, subsequent Coulomb blockade 
peaks are due to the filling of consecutive single-particle 
states 1]. Once a state is filled, it remains so. Al- 
though this picture successfully describes various trans- 
port properties of weakly-coupled meso- and nanoscopic 
systems, there has been much interest in identifying sit- 
uations in which the orthodox picture fails and a "dy- 
namical" behavior of the occupations of single-particle 
orbitals emerges. 

One of the earliest pertinent examples was pointed out 
by Kuznetsov et al. |2( who considered the filling of local- 
ized states in a barrier. They have shown that as the gate 
voltage increases, a localized state may first fill and then 
vacate, once a different localized state is occupied. This 
behavior is manifested in the conduction through the bar- 
rier by the reappearance of the same conduction peak. A 
new wave of interest was motivated by the correlations 
observed in the transmission phase through a quantum 
dot . An attractive explanation for these correlations is 
that a number of successive transmissionpeaks through 
the dot is carried by the same state 0, H, |(|. This, of 
course, requires the population of this state to be repeat- 
edly cycled. Various mechanisms, which lie beyond the 
orthodox picture, have been proposed 

In this paper, we show that repeated filling of a single- 
particle state is in fact a rather generic phenomenon. 
We find that it occurs whenever there exist traps in 
the system, either by accident or by specific design of 
a quantum-dot structure. This phenomenon has been 
seen by Lindemann et al. \J\ and in tailored structures 
by Johnson et al. @ and Kobayashi et al. |j| ■ Consistent 
with these experiments, we observe that repeated filling 
of a given single-particle state as function of gate volt- 
age can be reflected in the conductance in many different 
ways, ranging from essentially no signature to sawtooth- 
or dome-like structures to asymmetric Coulomb-blockade 
peaks. The underlying mechanism termed charge sensing 
in Ref. ||, is based on the capacitive coupling between 



the traps and the conducting channel. 

We first obtain our results within a reduced, exactly 
solvable model which captures the essential features of 
the phenomenon. As a by-product we show how the 
mechanism previously proposed by Silvestrov and Imry 
6] is a limiting case of our more general approach. Re- 
lying on numerical results obtained by a density-matrix 
renormalization-group (DMRG) method [13, EH- we sub- 
sequently discuss how our results are modified when re- 
laxing various restrictions of the exactly solvable model. 

Model. — Motivated by the experiment of Ref. 0, we 
consider a reduced model of the charge-sensing setup as 
shown in Fig. ^ A quantum dot is coupled to two leads. 
This connected quantum dot is coupled electrostatically 
to a disconnected dot in its vicinity. Within our reduced 
model, we make the following assumptions: (i) There 
is no tunneling between the disconnected dot and the 
leads or the connected dot. (ii) Transport through the 
connected dot is carried by a single state, (iii) The dis- 
connected dot may have many single-particle levels. We 
emphasize that this model may also represent a quan- 
tum point contact with nearby traps, as well as transport 
through a single quantum dot in which one state is much 
more strongly coupled to the leads than the others. 

Our reduced model is defined by the Hamiltonian 



H = H, 



clots 



Hi 



loads 



(1) 



Here, the Hamiltonian Hdots of the dots involves the level 
energy e of the connected dot (with creation operator ad) 
and e m of the disconnected dot (with creation operator 
d^j) together with the charging energies Ud and U m for 
the disconnected dot and for the mutual capacitive cou- 
pling of both dots. For simplicity, we will assume these 
charging energies to be equal. For spinless electrons this 
gives ffdots = ea^a+ Y, m Li e m d\ n d m + U[n d (n d - l)/2 + 
n a nd], where n a = a^a and n d = 53m=i ^m^m- The lead 
Hamiltonian is H ic£lds = Y,k,\=L,R E k\c\ x c k \ and the 
tunneling between lead and connected dot is described 
by -ffmix = Y.k\ tkxa^Ckx + h.c. 
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FIG. 1: Conductance (full line) through a single connected 
state with ten disconnected states as function of chemical po- 
tential, (a) Intermediate-coupling case: e = —1.08, t = 0.2, 
U = 0.16 and e m = -1.2 + 0.02m (m = 1, . . . , 10). (b) 
Strongly-connected case: e — — 1, t = 0.5, U = 0.08, 
e m = —1.24 + 0.04m. Dotted curves: Conductances for fixed 
occupations n d = 0, . . . , 10 of the disconnected states. As rid 
changes with fi, the conductance (full line) switches accord- 
ingly. Deviations in the tails stem from the finite band width. 
Inset: Schematic of the connected and disconnected dots. 



Generally, Eq. is a many-particle Hamiltonian. 
Note, however, that there is no charging term for the con- 
nected dot since it can be only singly occupied. When 
combined with the assumption of no tunneling to and 
from the disconnected dot, this fact makes our reduced 
model exactly solvable. Since the occupations (c?J„d m ) of 
the disconnected-dot states can only take on the values 
or 1, we can treat these operators as c- numbers in the 
Hamiltonian Eq. Q. Thus, the Hamiltonian can be re- 
garded as a set of 2 d single-particle Hamiltonians, one 
for each possible set {n m } of occupation numbers d\ n d m , 

H {n m } = (e + U n d )a^a + J~) e m n m 

m 

+ Un d {n d -l)/2 + H leaAs + H miK . (2) 

We first treat the zero temperature limit. The corre- 
sponding thermodynamic potentials at a given chemical 
potential y, are 



{n m } 



= H 



{n m } 



/z(cr a 



n d 



k,X=LM 



A\ c k\)- (3) 



The ground state occupation {n m } of the dots can now 
be found by determining the configuration with the low- 
est thermodynamic potential. Since 0/ n \ is a single- 
particle thermodynamic potential, this can be done by 
calculating Sl{ nm } for ah 2 Nd possible configurations. 

Analytical treatment. — Focusing on the essential 
physics, we present a full analytic treatment of the ther- 
modynamic potentials in Eq. (j3J for a single disconnected 
state of energy e±. The extension to several discon- 
nected states is straightforward. For a single discon- 
nected state, there are two different possibilities rid = 



and n d = 1 with corresponding thermodynamic poten- 
tials Slo an d As the chemical potential (gate volt- 
age) /x increases, the disconnected state will eventually 
be filled at /i = /i SW itch when the condition f2o — ^i is 
fulfilled. This switch in the occupation of the discon- 
nected state is accompanied by an abrupt change in the 
occupation of - and hence the conductance through - 
the connected dot. It is this general phenomenon that 
is referred to as charge sensing. We now proceed with a 
quantitative analysis. 

The Hamiltonians Hq and Hi, associated with 
empty and occupied disconnected state, are single- 
particle Hamiltonians with eigenstates \tpj) and 
and eigenenergies {e"} and {e}}- Thus, we can 
define the density of states on the connected dot 
z/0(!)( E ) = X\<^? (1) |a + c#? (1) }*(e - e° (1) ), and on the 



leads Af°«(0 =E,-<^° (1) |E 



„t 



We now express the relevant thermodynamic potentials 
as 



n nd = (ei-v)6 nd ,i+ / de(e- f i)(^(e)+Af nd (e)). 

J — oo 

Assuming that the density of states in the lead varies 
only on scales large compared to both U and the width T 
of the connected state, we can write v°{e) = (r/27r)((e — 
e) 2 + (r/2) 2 ) _1 and v 1 {e) =i/°(e-U). Moreover, a finite 
number of disconnected states will leave the continuum 
of lead states essentially unaffected so that Af° = A/ 1 . 
Thus, we have 

fii-fto = ei-fx+ r de(s-[i)(v (s-U)-v (s)) (4) 



and obtain 



Or — Oq 



d - fi + U/2 + 
H-U 



e 2U 
■ arctan 



arctan 



2(fi-U -eY 
f 



4n { ^-ef + {T/2f 



(5) 



upon performing the integration. 

We first consider the limit of a weakly- connected dot 
for which the distances of the connected dot energies 
from the chemical potential e — /i (for n d = 0), and 
e + U — ji (for n d — 1) are large compared to the 
level width T. Then, Eq. (J5J simplifies to ill — ilo — 
ei — e + (r/27r) ln(|/z — e — U\/\fi — e|), in agreement with 
the many-body perturbation theory result of Silvestrov 
and Imry 0. This implies that the switching occurs on 
the Coulomb-blockade plateau ((n a ) is an integer); hence 
there is no effect on the conductance at ^switch- However, 
the switch leads to repeated appearances of the same 
Coulomb-blockade peak [f| . 
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We now turn to the situation when the broadening T 
is comparable (intermediate coupling) or larger (strongly- 
connected) than the distances of the dot energies from the 
chemical potential /i. In these cases, the perturbation 
theory of Ref. 6] fails, while our general solution Eq. 
(J3J) still applies. Specifically as /i sweeps across Switch 
from below, the occupation of the connected-dot state 
decreases abruptly from n a = 1/2 + arctan[2(/! SW i tc h — 
e)/r]/7rto n a = 1/2 + arctan[2(/! SW i tc h — e — U)/T]/w. In 
the limit of a strongly-connected state, Eq. leads to 
the explicit solution ^/.switch — e± + U/2. 

Instead of the changes in the occupations n a , we focus 
directly on the experimentally more accessible (dimen- 
sionless) conductances g of the quantum-dot structure. 
For the specific case of a connected quantum dot with two 
symmetrically-coupled single-channel leads, the Friedel 
sum rule can be exploited in the usual way to derive the 
relation g = sin 2 (-7m Q ) 14]. Thus, the jump in the oc- 
cupation n a at /i SW itch translates directly into a jump in 
the conductance (unless An a = 1, which happens for a 
weakly-connected dot). 

Numerical results. — Representative traces of the con- 
ductance as a function of fi are shown in Fig. ^ These 
plots are based on a generalization of the above analytic 
results to the case of an arbitrary number N c i of dis- 
connected states. The most striking behavior occurs for 
intermediate coupling where the broadening T is compa- 
rable to or slightly smaller than the charging energy U. 
In this case shown in Fig.QJa), one observes the appear- 
ance of new asymmetric peaks in the conductance trace. 
In the absence of disconnected states, there would be 
only a single conductance peak due to the single level of 
the connected dot. In the presence of the disconnected 
states, the occupation of the connected state decreases 
abruptly whenever a disconnected state is filled up. Thus 
the sharp jump in the conductance is downward (upward) 
if it occurs on the rising (falling) side of the conductance 
peak of the connected level. This leads to the appear- 
ance of new trap-induced peaks in g whose asymmetry 
arises from the abrupt jumps. Indeed, hints of this be- 
havior have recently been seen in experiment and were 
attributed to charging of disconnected states Q. 

For strongly-connected dots (U <C T), the jumps in the 
occupation n a of the connected states which are associ- 
ated with charging of disconnected states are typically 
small compared to one. In this case, also the abrupt 
changes in the conductance are small compared to the 
conductance itself, leading to a characteristic sawtooth- 
behavior of g as function of chemical potential. This is 
shown in Fig. ^b). When the connected level is close 
to half filling, a typical dome shape is observed. This is 
very similar to the behavior seen in recent charge-sensing 
experiments on a quantum point contact monitoring the 
charge in a disconnected dot ||. 

In experiment, the gate voltage does not affect the con- 
nected and disconnected levels in the same way. Using 
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FIG. 2: Conductance as function of overall chemical potential 
fi and external gate voltage V g for the disconnected states (see 
inset in Fig. 0. Connected state: e = 0; disconnected states: 
e m = -0.55 + 0.05m (m = 1, . . . , 5), U = 0.1 and t = 0.2. 

different gates one can even manipulate the levels in- 
dependently |8(. Such experiments correspond to non- 
vertical trajectories in the (V g ,/j.) plane where y, is the 
overall chemical potential and V g is assumed to affect 
the disconnected states only. The resulting intricate pat- 
tern of the conductance in the (V g , fj,) plane is shown in 
Fig. El 

Temperature leads to very interesting behavior of the 
conductance, even when kT <C T. In the latter regime, 
the sharp jumps of the conductance are broadened by 
temperature, while the other side of the conductance 
peaks is smooth on the scale T and thus insensitive to T. 
This results in a very asymmetric temperature broaden- 
ing of the sawtooth- like peaks which was indeed observed 
in experiment Q. These considerations can be quanti- 
fied by noting that close to the charging of a disconnected 
state the system can be well approximated by a two-level 
system, with the two levels corresponding to — and 
nd = 1 (or, more generally, — N and n& = N + 1). 
The occupation of the connected state is then given by 

[§ + arctanffie-" + [\ + arctan 2{f *~ U) ]e~ Ul 



where uj nd — fl nd /kT. Since the entropy is governed 
by the lead states which are unaffected by the change 
in nd, one expects the entropy terms in fl to be equal, 
i.e., TSi = TSo, and thus to cancel out from Eq. 
Expanding Oi(/i switc h + Sfi) - tt (^switch + Sfi) ~ kT in 
<5/i, one obtains kT ~ <5/i(l — [arctan(2 (/Switch — U)/T) — 
arctan(2/i sw itch/r)]/7r. Therefore, the abrupt change in 
n a and hence g is smeared by temperature over a range 
kT in /i. The resulting T dependence of the conductance 
is depicted in Fig. EJa). 

A qualitatively similar effect occurs when the discon- 
nected state is broadened, e.g., by coupling it to an ex- 
ternal reservoir different from the current-carrying leads. 
This situation is no longer amenable to an exact solu- 
tion since the occupation of the "disconnected" state can 
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FIG. 3: (a) Temperature dependence of the conductance as 
function of chemical potential (connected state: e = —1.3; 
disconnected state: ei = —1.26; U = 0.4; t = 0.3). The 
different curves correspond to nkT — 0,0.1,0.2,0.3,0.4. (b) 
Occupation of the connected state n a (empty symbols) and 
"disconnected" state n,j (filled symbols) for several values of 
couplings as the "disconnected" state opens up additional 
non-current carrying lead. Lines: results of the exact solu- 
tion of Eq. 10; symbols: DMRG calculations for different 
couplings of the "disconnected" state (circles t = 0, squares 
t — 0.1, diamonds t = 0.2). (c) Conductance of connected 
state in (b). 



now differ from or 1, and we resort to a numerical 
DMRG method. (We have checked that the DMRG re- 
produces our exact solution for vanishing broadening, see 
Fig- ED 5 ).) We find that the main effect of a finitely cou- 
pled "disconnected" state is to smear the sudden jump in 
the occupation of the latter, similar to the effect of finite 
temperature. Correspondingly, the effect of the broad- 
ening on g is very similar to that of temperature, as can 
indeed be verified by comparing Figs.^a) and[3fc). 

This observation enables us to qualitatively under- 
stand the physics of two connected states of the same 
dot with T± ^> IV Indeed, the level occupations are in- 
sensitive to whether the levels are coupled to the same or 
different leads. Fig. O^b) then applies and is consistent 
with very recent works on such setups emp loying Hartree- 
Fock and numerical RG methods |l2|, [l3| . Evidently, the 
way these states are coupled to the leads is crucial for 
the conductance, due to interference effects such as Fano 
resonances 0,0. 

The insight gained from our exact analysis involving 
a single connected state may be used to make quanti- 
tative predictions concerning several connected and dis- 
connected states, including charging energies for the 
connected-dot states. Although it is then difficult to 
calculate the exact switching point for a particular dis- 
connected state, one nevertheless predicts that for chem- 
ical potentials /i (or V g ) immediately before or after 
the switching, n a j{n ± 0) = n a j(fi T- U) and therefore 
g(fi ± 0) = g(fi T U). (Here, n a j denotes the occupation 
of the j'th connected state.) This behavior is illustrated 
in Fig. 2| for two connected and two disconnected states, 




FIG. 4: Occupation for two connected states (t = 0.5) at 
ei = —1.1 (gray circles) and £2 = — 1 (gray squares) and two 
disconnected states ei = —1.12 (line), = —1.02 (dotted 
line) as function of fi calculated using DMRG for U = 0.3. To 
demonstrate that right after the switch n a (/x) = n a ([i — U), 
we plot n a (fi + U) for values of fx prior to the switch (i.e., 
-1.2 < < -1, and -0.48 < ^ < -0.28) shown by white 
symbols. A clear overlap between the white and gray symbols 
is seen. 



based on a DMRG calculation. This prediction should 
be very useful for analyzing experimental data. For any 
abrupt jump due to charging of a disconnected state, the 
conductance satisfies this relation. By contrast, if the 
jump in the conductance is due e.g. to noise no such cor- 
relation is expected. 

In conclusion, it is interesting to speculate that the 
charge-sensing physics discussed theoretically in this 
paper may occur generically in relatively well-coupled 
chaotic quantum dots. By the nature of the Porter- 
Thomas-distribution of lead-induced level broadenings, 
there will be a significant number of narrow levels in ad- 
dition to broader levels, in particular for systems with 
time-reversal symmetry. Charge sensing and patterns 
not unlike Fig.^b) may thus be important ingredients in 
explaining the large-scale structure of Coulomb-blockade 
sequences observed in such systems [l5j |. 
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